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ABSTRACT 

The atomic oxygen concentration in the altitude range 130-240 km 
has been determined through the use of a quadrupole spectrometer with 
a strongly focussing ion source. The instrument is used in the Thermo- 
sphere Probe in a manner that greatly increases the proportion of mea- 
sured oxygen ions that have not experienced a surface collision and 
permits quantitative evaluation of surface recombination and thermali- 
zation effects which inevitably enter all spectrometer determinations. 
The data obtained strengthen the concept that consideration of surface 
effects is significant in quantifying spectrometer measurements of re- 
active gases, and tend to be in agreement with von Zahn’s recent 
results. 
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A THERMOSPHERE COMPOSITION MEASUREMENT USING A 


QUADRUPOLE MASS SPECTROMETER WITH 
A SIDE ENERGY FOCUSSING QUASI-OPEN ION SOURCE 

1. INTRODUCTION 

The measurement of atomic oxygen by means of rockets has posed difficult 
instrumental problems because of the relatively '.indefinable but dominant gas/ 
surface interactions experienced. At satellite altitudes where 0 greatly pre- 
dominates over 0 2 , the 0 concentration can be measured relatively easily by re- 
combination of the 0 into 0 2 in the mass spectrometer ion source chamber. 
However, in the lower thermosphere, where 0 and 0 2 coexist in delicate equi- 
librium, this approach is inadequate. We believe however that a satisfactory 
approach is one that employs a mass spectrometer that ionizes and measures a 
significant number of 0 atoms that have not experienced significant surface col- 
lisions in the ionization chamber. As a result, in the data analysis, a model can 
be used, which takes into account both particles which have and have not experi- 
enced surface interactions. In the experiment discussed in this paper, an instru- 
ment which was equipped with a quasi-open source and provided adequate side- 
energy focussing, and thus met the required conditions was employed. Because 
of the focussing, we were able to consider measurements at both near-zero and 
90° angles of attack of the ion source. The system thus permitted evaluation of 
the relative gas/surface effects and made possible, we believe, unambiguous 
measurements of the atomic oxygen concentration in the thermosphere. We 
present the vertical profiles of the partial densities of 0, 0 2 , and N 2 obtained 
over the altitude region from 150 to 280 km. In addition, the techniques employed 
to analyze the data using a time-independent absorption model are discussed. 
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2. INSTRUMENTATION 


In this experiment, a Thermosphere Probe (TP) instrumented with two mass 
spectrometers, a quadrupole and an omegatron, was launched from Wallops Is- 
land, Virginia, on August 26, 1966, to study the earth's atmospheric parameters 
in the altitude region between 120 to 350 km. The launch vehicle was a Nike- 
Tomahawk. It was launched at 14:31 EDT as one of a series of five launches 
scheduled to measure the extremes of the diurnal variation of atmospheric tem- 
perature and density. 

The Thermosphere Probe (TP) which is launched in a clam-shell nose 
cone (Figure 1) and then ejected into the lower thermosphere with an end over 
end motion contained the quadrupole and omegatron mass spectrometers from 
which the data discussed herein were obtained. The spectrometers, mounted in 
the TP wdth their ion sources centered at each end, were kept in vacuum until 
an altitude of about 150 km was reached, where the automatic breakoff devices 
exposed the instruments to the ambient atmosphere. A short time later, mea- 
surements were initiated by turning on the filaments. The omegatron spectrom- 
eter was located behind an orifice in the conventional manner reported by 
Spencer, et al., 1965. The quadrupole mass spectrometer, however, utilized a 
quasi-open side energy focussing ion source to reduce the effect of gas-surface 
interactions on the atomic oxygen measurements. The source used was similar 
to that employed with the Explorer 17 mass spectrometers (Spencer and Reber, 
1964; Reber and Hall, 1966) but its electrodes were made of gold plated stainless 
steel and the electrode potentials were adjusted for rocket flight conditions so 
that free streaming neutral particles could be analyzed in the mass spectrometer 
without energy and momentum discrimination. 
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A detailed view of the ion source is given in Figure 2. Electrons are 
emitted from a hot tungsten-rhenium filament and are accelerated by an electric 
field generated by a potential of about 75V applied between the top surface and 
the filament. A permanent magnetic field focuses the electrons, forcing them 
to move along an arc as indicated in the figure. Neutral gas particles are ion- 
ized by electron impact along the electron path. A positively-biased highly 
transparent (95 percent) repeller above the top surface, and an ion lens system 
below the center opening of the top surface provide an electric field which ex- 
tracts ionized particles from the primary ionization region indicated in Fig- 
ure 2. The lens system produces a magnified image of this region on the source 
base. A portion of this image is projected through the entrance aperture (marked 
in Figure 2 with an X) into the quadrupole analyzer section. This ion lens geom- 
etry is designed, at the expense of sensitivity, to assure that the transmission 
of ions resulting from ambient neutrals arriving in the source region is inde- 
pendent of the initial velocity and incident angle of the neutral particles entering 
the source region for rocket velocities up to about 2 km/sec. In this manner, 
neutral particles entering the ionization region directly from the stream are 
detected thus fulfilling the requirement that particles not experiencing a wall 
collision will be measured. Because some particles that have reacted with the 
walls and have been thermalized also inevitably are detected the ion source is 
referred to as quasi-open. The ratio of the concentration of free stream ions 
to thermalized ions is a strong function of the angle of attack and provides the 
basis for the primary analysis considerations of this paper. 

The quadrupole analyzer section was operated in a 100% ion transinission 
mode which resulted in flat-topped mass peaks and permitted acceptable 
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tolerance for normal drifts of the operating voltages and frequencies. Selection 
of masses 40, 32, 28, and 16 was accomplished by varying the rf-oscillator fre- 
quency and keeping the rf- and dc-voltages constant. Two other modes of op- 
eration included first, a total density measurement for which the dc-voltages 
were disconnected from the quadrupole rods and the rf-frequency was tuned to 
the mass 16 value. Transmission of all ions with atomic masses above, ap- 
proximately, mass 13 was thus assured. The other mode was a zero-transmission 
mode for ion detector calibration and noise evaluation in which the ion flow was 
interrupted by changing the ion lens voltage and tuning the rf-oscillator to a 
frequency corresponding to mass 25. The ions were detected at the end of the 
analyzer section with the aid of a secondary electron multiplier. A multiplier 
gain of 10 s provided adequate current sensitivity. The anode current of the 
multiplier, varying from 10' 12 to 10" 7 amperes was converted to a voltage 
range of 0 to 10 volts by a logarithmic electrometer amplifier. 

During ejection from the nose cane the TP was caused to tumble end-over- 
end with a period of 2.62 sec. The quadrupole was tuned sequentially to the 
masses 32, 16, 28, 40, and zero for 2.8 sec each, corresponding to about 1.05 
tumble periods, which permitted measurement through a complete tumble period 
for each species. During the zero ion current mode, the electrometer was cali- 
brated at three current levels: 10' 6 A, 10~ 9 A and zero current. Near apogee 
the sequence was lengthened to include the total current measurement. 

Before flight the quadrupole spectrometer was pressure-calibrated simul- 
taneously with four omegatrons against three Bayard-Alpert (BA) ionization 
gauges which had been calibrated against a McLeod gauge. The calibration with 
the BA gauges gave meaningful results for nitrogen and argon, but since BA 
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gauges generate carbon oxides and exhibit a memory effect when operated in an 
oxygen atmosphere, only a relative oxygen calibration between the omegatron 
and the quadrupole was obtained. Thus, during the oxygen calibration, the BA, 
gauges were turned off. The contribution of dissociatively ionized oxygen to the 
atomic oxygen measurement was determined from the mass 16 to mass 32 cur- 
rent ratios. For absolute oxygen sensitivities, the ionization cross section ratios 
cr( 0 2 )/cr ( N 9 ) and cr (0) /o-(N 2 ) (Kieffer and Dunn, 1966, and Fite and Brackmann, 
1959) were used together with the experimental nitrogen sensitivity value. 

The quadrupole system had two in-flight difficulties: (a) post-launch com- 

parison with the omegatron made it apparent that the sensitivity for N 2 had de- 
creased by approximately a factor of four, and (b) the rf-oscillator failed to 
tune to the proper frequency for mass 40. Subsequent studies in the laboratory 
on likely causes for the sensitivity change were { 1) a change in the bias voltage 
of one of the ion source lenses which increases the magnification of the lens and 
thus decreases the sensitivity of the instrument, or (2) a change in the electron 
multiplier gain. The most probable cause of oscillator detuning at mass 40 was 
a failure of the reed relay used to change the oscillator frequency by capacitance 
switching. 

In the subsequent analysis, it was assumed that ion source focusing was not 
adversely affected by the change in sensitivity and that the omegatron nitrogen 
measurement provided an adequate basis for establishing the level of in-flight 
quadrupole sensitivity. 
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3. ANALYSIS OF DATA 


The tumbling motion imparted to the TP produces a periodic enhancement 
and rarefaction of the particle density in the ionization region of the quadrupole 
ion source. This is predicted qualitatively by the thermal transpiration equation 
for a moving gas (Schultz, et al., 1948) , which applies to all closed and idealized 
ion sources. For the particular geometry of the quadrupole ion source, a quanti- 
tative analysis made by Niemann and Kreick (1966) applies. 

During this flight the roll rate (period: 1.075 sec) and the tumble rate 
(period: 2.62 sec) of the TP resulted in a coning angle of 4.1 degrees. The sun- 
aspect sensor data and trajectory data permitted a continuous measurement of 
these angles, and thus the angle between the ion source normal and the velocity 
vector. Data were taken for 225 sec of flight time yielding 22 sets of data points 
for each mass. Since each mass was observed for 1.05 tumble periods, data 
were obtained at all available angles of attack. Figure 4a shows the envelope of 
the minimum and maximum angles experienced throughout the flight, and 4b 
shows the variation of instantaneous angles of attack with time for a short period. 

Using in-flight omegatron data and published ionization cross sections, 
particle densities in the ionization region were calculated. The cross section 
ratios used in this analysis were: cr (0)/cr(N 2 ) = 0.786 (Fite and Brackmann, 
1959); cr (o 2 ) /cr (N 2 ) = 1.0 (Kieffer and Dunn, 1966). Plots of source densities 
for the three gases (molecular nitrogen, molecular oxygen and atomic oxygen) 
versus time are shown in Figures 5,6, and 7. In each figure, the same densities 
are shown for three distinct orientations: curve (a) minimum angle of attack 
attained during a given tumble, curve (b) 90 deg angle of attack, and curve (c) 
maximum angle of attack. The atomic oxygen densities were also corrected for 
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the contribution of dissociately and doubly ionized molecular oxygen which, 
because of the relatively high energy of the ionizing electrons and the strong 
focusing properties of the ion source, amounted to 50 percent of the indicated 
0 2 density. However, because of the rapid decrease of the ambient molecular 
to atomic oxygen ratio with increasing altitude, the contribution is very small 
over much of the higher altitude region. 

To relate the measured particle density in the ionization region to the 
ambient particle densities, it is necessary to account rigorously for the gas 
dynamic characteristics of the ion source. Although the ion source is relatively 
"open", there are significant contributions to the density in the ionization region 
at some angles of attack by particles which have had surface collisions and, 
because of the focussing properties of the ion source, by those which have not. 

In detailed consideration, four groups of particles, each group having a 
distinct history, constitute the instantaneous particle population of the ionization 
region: (a) particles which stream directly through the ionisation region without 

having made any collision with the instrument surfaces, (b) particles which enter 
the ionization region after one bounce with the instrument surface, (c) particles 
which make multiple surface collisions before entering the ionizing region, and 
(d) particles which are slowly emitted from the surface of the instrument as a 
nearly constant background contribution. The ambient particles which have 
undergone surface collisions (single and multiple bounce) are further subdivided 
according to the nature of the surface collision; i.e., specular reflection, diffuse 
reflection, degree of thermalization, degree of absorption; the relative contri- 
bution of these various classes of particles are evaluated as a function of the 
variable speed ratio and angle of attack. Specifically, for the atomic oxygen 
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data reduction, a gas-surface interaction model is chosen to fit the data at the 
minimum angle of attack and at 90 deg angle of attack. Assuming that the ab- 
sorption and recombination properties of the surfaces change slowly in time 
compared to the transition times of the gaseous particles, relations between the 
ambient particle density and the velocity and angle of attack dependent particle 
density in the ionizing region are obtained by integrating the Maxwell-Boltzmann 
velocity distribution function in velocity space over the proper limits for the 
geometry shown in Figure 8. The final expressions are given in Equations 
(3.1) through (3.5) , which follow. The definitions of the terms used in these 
equations are: 

n 0 = ambient particle number density, 

n A = Number density of neutral particles in primary ionization region, 

Tj = instrument temperature, 

T 0 = ambient neutral particle temperature, 

K = permeability of the repeller screen, 

8 , = acute angle between the outer edge of the deck and the deck normal 

measured from the ion source, 

8 2 ~ acute angle between the edge of the orifice and the deck normal 
measured from the ion source, 

S = speed ratio, I u|/(2kT 0 /m)^ 
where u = vehicle velocity, 

k = Boltzmann's constant, 1.38 x 10 -23 J/°K, 
m = molecular weight of the specie measured, 

(2kT 0 /m) ,/j = most probable speed in a gas at temperature T Q . 

S z = S cos a, 
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where a. — angle of attack between u and the deck normal. 

3 r = S sin a, 

n d j r = contribution of particles which stream through the ionization 

region, without collisions with the top surface prior to ionization. 
n spec = contribution of particles which stream through the ionization 

region after a single collision, with the top surface if all particles 
are specularly reflected, 

n diff ~ contribution of particles which stream through the ionization 

region after a single collision with the top surface if all particles 
are diffusively reflected (according to the cosine law) without 
thermal accommodation, 


n dif f , acc = contribution of particles which stream through the ionization 

region after a single collision with, the top surface if all particles 
are diffusively reflected with thermal accommodation to the tem- 
perature of the surface, 

n orif = contribution of particles which stream through the ionization 

region after many collisions with the surfaces behind the orifice. 
Complete thermal accommodation and an ideal cavity are assumed. 


dir 


K 


9 [l +■ erf (3 Z ) + exp(-S z 2 ) M 0 (S z , l9 1 ) + H(S z ,S r ,0 1 )] 


(3.1) 


where 
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M 0 - cos 6 j exp (S 2 cos 2 ^) [1 - erf (S z cos0 1 )] , 


9 



M l = (2 + S z 2 COs2 ^i) M 0 - S z COs2 ^! ^ - 

Mi = f {[S z 2 cos^j ,+ |(4i - 3)] M i . 1 - ~(2i - 3) M^}, i > 2. 

Spec Kf nr 

— = - |exp(-S2) [M 0 (-S z ^ 2 )-M 0 (-S z ^ 1 )]fH(-S z> S r ^ 2 )-H(-S 2 ,S r ^ l) - . 

(3.2) 

— = | [l + erf (S z )] [cos (9 2 - cos 6 x ] . (3.3) 


where 


n d iff , acc 


n 


0 



[cos 6 2 


cos 6 j] , 


f(S z ) - exp(~ S z 2 ) + fff S z [l + erf(S z >] 

n o r i f K ^0 

~ = 2 T f (SJ (1 “ cos9 2 ) . 

0 l 


(3.4) 


(3.5) 


If we let F x be the probability that, a particle striking the top (deck) of the 
ion source will be absorbed, and F 2 be the probability that a particle entering 
the ion lens region through the hole in the deck will be absorbed, we can write 
an equation which relates the source particle density to the ambient density in 
terms of the surface parameters and the gas kinetic relation given above. 


n. d . 

n - = n o 1~7T~ r + O ~ F i) 

' 0 


spec 

a, — — + a 
L 1 "O 


‘di f f , acc 


^ n. 3 n 


(3.6) 


(1 - F a ) 


o r i f 


n. 


+ n. 


bkg 
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The last term, n bkg , is added to account for the contribution of background 
gases which are absorbed on the ion source surfaces during ground calibration 
and preparations, and are then emitted during the flight. Of the particles col- 
liding with the deck, a. 1 is the fraction which shows specular reflection, a 2 dif- 
fusive reflection without thermal accommodation, and a 3 diffusive reflection 
with thermal accommodation. These fractions depend on the species and are, 
therefore, different for 0, 0 o , and N 2 . 

The assumption that the particles reflected from the deck can be classified 
into three distinct groups is a simplification. More realistically, particles come 
off the deck according to a distribution which is somewhere between specular 
and diffusive with a characteristic temperature lying between their original and 
that of the deck. For a better approximation, an additional expression Is pre- 
sented which combines the equations for the diffusively reflected., thermally ac- 
commodated and non-thermally accommodated particles into one equation by 
using a modified velocity distribution function for the reflected particles. In the 
derivation, it is assumed that the incoming particles lose an amount of energy 
to the deck which is proportional to their initial energy. The expression which 
replaces the surface term 


spec 


Miff 


+ a. 


-i- a~ 


d i f f , a c c 


is 


n 


ref 


i (O 


K 

2 f(S') 


T„ 


(cos Q , 


cos 6 j ) 


1+ e{~ - l) [1 + erf (S;)] , (3.7) 


where 
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1 - e 


Z 



The term e is the accommodation efficiency defined: 


V-U 


hi 2 


c o 2 

I -* — 1 1 2 0 I — * | 2 

V-U I | V 

c 2 


where 

u - probe velocity, 

v = particle velocity in the stationary coordinate system, 
v' = final particle velocity after collision, 

C s = i/2kT s /m, most probable speed in a gas at the temperature of the 
surface, 

C Q - |/2kT 0 /m, most probable speed in a gas at ambient temperature. 

A derivation of Equation (3.7) is given in Appendix A. In the subsequent reduction 
of the data both models were compared. 


3.1 Molecular Nitrogen 

The ambient nitrogen density .may be computed by using Equation (3.6) for 
any orientation provided that TP velocity, TP angle of attack, accommodation 
coefficients and ambient temperatures are known. The constant background 
contribution to the source density is eliminated by taking the differences between 
the output at some known angle of attack and the minimum measured output as 
in the omegatron experiment data analysis. Equation (3.8) takes this procedure 
into account and follows from Equation (3.6). Here it is assumed that the surface 
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reflects all particles diffusively (a x = 0), that no adsorption occurs (F t - 0, 
F ? - 0), and that the accommodation coefficients are independent of angle of 
attack. 


An i = ^o( N 2>1 A 


dir 


+ a„ A 


‘di f f 


n 0^2) 2 n o(^ 2 ) 


i- (1-a ,) A 


l diff,acc 

n Q (N 2 ) 


+ A 


li • 

o r i t 
n 0( N 2) 


(3.8) 


where 

An. - n. (a) - n. (a ), etc. 

a = chosen angle of attack 
a = maximum angle of attack. 

rn ax 0 

The ambient nitrogen density may also be computed from Equation (3.9) in 
which Equation (3.7) is used to estimate surface contributions to the ion source 
density. 


An. 


n 0 (N 2 ) iA 


re f 1 “onf 

+ A TTT" (e) + A 


n o( N 2 ) n o( N 2 ) 


n o(^ 2 ) 


(3.9) 


The computed source density-to-ambient density ratio depends on the 
ambient gas temperature which during this experiment was measured by the 
omegatron. The uncertainty in the thermal accommodation parameter can be 
reduced by applying Equation (3.8) at different angles of attack since the angular 
dependence of each term in Equation (3.8) is different. A suitable second de- 
termination of the accommodation coefficient is obtained when the angle of attack 
is 90 degrees, where the direct-streaming contribution is relative^ large, and 
all the terms in Equation (3.8) or (3.9) depend only weakly on the ambient tem- 
perature. Again, the differences between these readings and the readings at the 
maximum angle of attack were used to eliminate the background. 
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Equation (3.8) was used to solve for the ambient nitrogen density, n 0 , by 
using the given parameters and the measured source densities from both upleg 
and downleg data, using source densities at different angles of attack. Unknown 
parameters, namely, accommodation coefficients and ambient temperatures, 
were determined by requiring that all calculations be consistent with a unique 
ambient density profile. To determine the accommodation coefficients of nitro- 
gen on the quadrupole ion source deck, the following approach was used. 

An ambient profile was calculated for each of two extreme assumptions of 
(a) complete thermal accommodation of the particles to the deck, and (b) zero 
thermal accommodation. Calculations for both minimum angle of attack and 
90 deg angle of attack were made. In this manner, four upleg density profiles 
for differing conditions were obtained, which are plotted in Figure 9. The same 
procedure was followed to obtain density profiles using the downleg data. 

A single density profile estimate was then determined from the four cases 
by assuming accommodation coefficients which were (a) constant for the entire 
measurement period, and (b) resulted in closest agreement between the available 
calculated profiles. The resulting single profile is shown in Figure 10, and 
corresponds to coefficients a, - 0.56 and a 3 = 0.44. By using the accommodation 
efficiency correlation with e =0.52, employed in Equation (3.7), a separate set 
of data were computed which are also shown in Figure 10. These data differ 
little from those resulting from use of the coefficients a 2 and a 3 . 

The downleg data obtained at 90 deg angle of attack yielded ambient densities 
in exact agreement with the upleg results. However, the downleg data obtained 
from measurements at minimum angle of attack yielded ambient density values 
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approximately 20% higher than the upleg values. It was concluded that a linear 
scale adjustment to force these data points to coincide with the rest was appro- 
priate, as will be discussed later in the paper. 


3.2 Molecular Oxygen 

The molecular oxygen data reduction procedure is similar to that followed 
for nitrogen. The computations are more complicated, however, because mo- 
lecular oxygen exists also on the surface of the instrument due to recombination 
of atmospheric atomic oxygen. To separate the contribution of the atmospheric 
molecular oxygen and the recombined atomic oxygen, atomic recombination must 
be taken into consideration. The expression for the total source (oxygen) density 
is: 


An i(° 2 ) = n o(° 2 ) A 


Mi i 


+ a' A 


‘di f f 


n o(° 2 ) 2 n o(° 2 ) 


+ ( 1~ a' 2 ) A' 


a i f f , ac c 

n o(°2) 


i- A 


o r i f 

v»7) 


(3.10) 


n 0 (°) 


n . 


Miff n diff 

A tttt + ( 1- a!) ) A rr? 
n 0 (°) v ^ n Q (0) 


1 “ R o r i f 
+ A 


n o(°) 


where -y denotes the probability of recombination after one collision and gener- 
ally depends upon particle flux, and R is the fraction of recombined atomic oxy- 
gen exiting the ion lens region given by: 

1 

R = — 

1 + p ( i - y) 

p is the probability of escape from the ion lens region after a collision, which 
is assumed to be the same after each collision (Niemann and Kreick, 1966). The 


15 



primes on accommodation coefficients denote the different degrees of accommo- 
dation involved for atomic and molecular oxygen. 

From these expressions, density profiles were computed for the extreme 
cases of complete thermal accommodation and of zero thermal accommodation. 
Several values of y and R were tried, but only those profiles resulting from the 
assumption of zero recombination of atomic oxygen had the proper altitude de- 
pendence for diffusive equilibrium. At higher altitudes, even a very low re- 
combination rate would have been apparent in the molecular oxygen measure- 
ment because of the high 0/0, ratio. It was concluded therefore that during the 
time of measurement, the atomic oxygen was either absorbed, reflected, or 
possibly formed another constituent by reacting with some contaminant which 
was adsorbed on the surfaces. The upleg results of the ambient density compu- 
tation similarly without atomic oxygen recombination are shown in Figure 11. 

The choice of the coefficients a 2 = 0.21 and a 3 = 0.79 for both orientations 
resulted in a single profile as shown in Figure 12, where downleg results are 
also plotted. The downleg data obtained from minimum angle of attack measure- 
ments also showed values which were 20% higher than the data obtained for the 
upleg and 90“ angle of attack measurements. A scale adjustment was made for 
these data points on Figure 12. 

3.3 Atomic Oxygen 

The analysis of atomic oxygen data is more complex because of the stronger 
chemical reactions it experiences with the ion source surfaces. Oxygen atoms 
can react in several ways with the instrument surface, the most probable re- 
actions appearing to be atomic recombination, permanent adsorption, and the 
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formation of gaseous constituents by reaction with previously adsorbed gases. 
Laboratory simulations of flight experiments and direct flight tests with omega- 
tron mass spectrometers have shown that atomic oxygen is strongly adsorbed 
on "clean 11 instrument surfaces for time periods which greatly exceed the mea- 
surement duration of this experiment. In all omegatron flight experiments of this 
kind, no oxygen was detected in atomic form for the sensors employed, that were 
constructed from either stainless steel or gold-plated stainless steel. 

To analyze the atomic oxygen data measured by the quadrupole it is, there- 
fore, necessary to consider adsorption phenomena. Since the specific surface 
characteristics are not well known, some simple assumptions were made to 
compute the extreme values of the possible ambient atomic oxygen density pro- 
files. As in the analysis of the 0 2 and N 2 data, the criterion that the assumptions 
made be consistent at all angles of attack again provides a powerful tool. 

It was established as discussed above that in the molecular oxygen data 
analysis that no significant atomic oxygen recombination occurred which con- 
tributed to the molecular oxygen measurement. In the atomic oxygen analysis, 
therefore, it follows that there is no loss of atomic oxygen due to recombination. 
By using the same procedure as in the reduction of the nitrogen data the equation 
for the source density difference is obtained. 


dir 
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where b 2 , b 3 , and b 4 are numbers equal to or smaller than unity. These numbers 
allow for fractional adsorption and the loss of atomic oxygen due to reactions on 
the surfaces to be included in the gas -dynamics model. Based on the results of 
the omegatron experiments, it is assumed that the particle contribution of the 
atomic oxygen streaming back out of the orifice to the ionizing region was negli- 
gible and that only direct streaming particles, and particles reflected from the 
top deck contributed to the atomic oxygen source density. This assumption re- 
quired that b 4 - 1, so Equation (3.11) becomes 


VO) = n Q (0) A 


n o(0) 


+ a_(l-b 2 ) A 


“diff 

n o(°) 


(1-a,) (l-b 3 ) A 


Miff , a c c 

VO) 


(3.12) 


The upleg ambient atomic oxygen density was computed using this equation 
for three extreme cases: (a) total adsorption on the deck (b 2 = b 3 = 1), (b) no 
adsorption or recombination on the deck with all particles reflected diffusively 
(zero accommodation), and (c) no adsorption or recombination on the deck with 
all particles reflected diffusively after complete thermal accommodation. The 
resulting six density profiles for minimum angle of attack and for the 90 deg 
angle of attack are shown in Figures 13 and 14. Again, to bring the minimum 
angle of attack results and the 90-deg angle of attack results into agreement, the 
coefficients were adjusted. 

Considering these results in contrast to the cases of nitrogen and molecular 
oxygen, the slopes of the atomic oxygen profiles obtained from the 90 deg angle 
of attack were lower than the slopes of the minimum angle of attack profiles. 
Solving for the accommodation coefficients (assuming they were constant) yielded 
close, but still unacceptable, agreement for the cases of no adsorption and full 
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thermal accommodation. This profile was essentially identical to profile C in 
Figure 13. The slope of this curve, however, was considerably less than that 
required by assuming thermal and diffusive equilibrium of the atmospheric 
atomic oxygen. 

As noted before, in a closed source where many wall collisions occur, 
atomic oxygen is completely adsorbed by gold plated stainless steel walls during 
the period of a rocket flight. However, the probability of adsorption after one 
bounce is unknown. Assuming that the ambient atomic oxygen was in thermal 
and diffusive equilibrium, a density profile was computed with an arbitrary mag- 
nitude but with a slope which corresponded to the ambient temperature obtained 
from the accompanying omegatron experiment. This theoretical profile was 
compared to the profile obtained from the measured data and the disagreement 
between them was attributed to surface sorption and possibly atomic recombina- 
tion forming compounds other than molecular oxygen. To be physically meaning- 
ful, the theoretical density profile so determined must lie between the curves 
A and C of Figures 13 and 14. The spread between the curves A and C in Fig- 
ure 14 is narrower than that in Figure 13, thus limiting the possible choice to 
that indicated by the dashed lines in Figure 14. The upper curve would be valid 
if total adsorption occurred at the low altitude, and the lower curve would be 
valid if no adsorption occurred at the high altitude extreme. The correct profile 
must then lie between the lines whose ordinates differ by about 16 percent. 

A theoretical adsorption curve was fitted to the curves in Figures 13 and 
14, and an improved ambient density profile was computed with the aid of this 
model. The simplest profile which could be fitted to the curves was a multilayer 
adsorption model derived by letting the adsorption characteristics of the first 
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layer of atoms on the surface be different from the characteristics of the higher 
order layers, and by using a steady state approximation. The derivation of this 
model is given in Appendix B and is similar to the Brunauer-Emmett-Teller 
adsorption Model (Flood, 1997). Only the final numerical expressions are given 
here. For the 90 deg angle of attack, 

r? 

i ' 90 ° P 2 - (.528 x 10 13 )F 1 + 6.89 x 10 26 ’ (--13) 

and for minimum angle of attack, 


r i 2 

F - .517 , (3.14) 

™ in B 2 + (3.23 x 10 ! - 3 )rj + 2.00 x 10 25 

where Fj is the incoming flux density of atomic oxygen (part/ cm 2 sec), and F is 
the fraction of incoming flux that is adsorbed by the deck. The coefficients a 2 
and a„ were arbitrarily chosen to be 0 and 1 respectively for this component. 

An upleg ambient density profile including adsorption was computed and is shown 
in Figure 15. It should be noted that this procedure is not unique and that dif- 
ferent surface adsorption models may yield similar atomic oxygen density pro- 
files. Therefore, no conclusion about the surface reactions was drawn from this 
lysis . 


In contrast to the upleg density profile the profiles computed from the down- 
leg measurement had the same slope as would be expected from a gas in diffusive 
and thermal equilibrium. A final downleg profile was, therefore, computed by 
solving for the coefficient a 2 and a 3 and assuming no surface adsorption. Agree- 
ment with the upleg data was obtained by choosing the coefficients for the 90° 
angle of attack a 2 - 0.34 and a 3 - 0-66. The downleg shift in apparent sensitivity, 
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observed for the molecule species near minimum angle of attack was a.lso found 
in this measurement. 

3.4 Temperature 

The kinetic temperature was determined by integrating the differential 
equation formed by the hydrostatic equation and the equation of the state of the 
idea .1 gas. 

f ° P( h '> g( h ‘) dh ' ,, , 

P( h o) 

T(h) ' Rp(h) + p (h) ^ h o) (3,1 ^ 

where 

g = gravitational acceleration 
/ 3 (h) ~ mass density at the particular altitude h 
h - altitude 

h Q = reference altitude where h 0 > h 
T(h) = neutral particle temperature at the particular altitude h 
R = universal gas constant (8.314 J/mole-deg). 

The technique described above is essentially the same as that used to com- 
pute the temperature profile from the nitrogen density profile measured with the 
omega/cron. All quantities in Equation (3.15) are taken from the quadrupole re- 
sults except T(h 0 ), which is taken from the omegatron data (Spencer et al., 1970). 
The resulting profile is shown in Figure 18. As expected, there is agreement 
with the nitrogen temperature profile of the omegatron. Molecular oxygen and 
molecular nitrogen profiles measured by the quadrupole are also consistent with 
an atmosphere in thermal and diffusive equilibrium, and the atomic oxygen density 
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profile has been constrained to yield the same temperature as that of the 
molecular nitrogen profile. 
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4. DISCUSSION 


Before conclusions can be drawn from these data, it is necessary to discuss 
the errors which contribute to the uncertainty in the measurements. The errors 
common, to the omegatron and the quadrupole experiments are those due to er- 
rors in trajectory, velocity, angle of attack, and calibration. They have been 
discussed by Spencer, et al. (1965), and for the particular case under considera- 
tion are expected to be as follows: maximum error due to position uncertainties 
is ±5 percent. The velocity uncertainty is smaller than ±1 percent, which results 
in a density error of about ±1 percent. The angle of attack, obtained by an opti- 
cal technique using the method described by Carter (1968), is believed, to be 
within ±5 deg, which results in a maximum error of ±10 percent at 90 deg angle 
of attack and less than ±2 percent at minimum angle of attack. The telemetry 
error is less than ±2 percent. 

By far the largest uncertainty is that of the pressure calibration. Both the 
quadrupole and the omegatron were calibrated simultaneously on a dynamic 
calibration system against several Bayard-Alpert (BA) ionization gauges, used 
as secondary standards. The absolute accuracy was believed at the time to be 
approximately ±25 percent for molecular nitrogen. However, present indications, 
and new knowledge about MacLeod gauge trap pumping effects would lead one to 
increase the calibration constants by constant factor of 1.2. Absolute calibra- 
tions for molecular oxygen could not be obtained because of the memory effect 
of the BA gauges in an oxygen atmosphere. Also, because an adequate technique 
does not exist, atomic oxygen calibration could not be performed. Relative cali- 
brations were carried out, however, for molecular oxygen and were automatic 
for molecular nitrogen. The relative calibration accuracy between the omegatron 
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and the quadrupole was considered to be within i5 percent but, because of the 
change in the sensitivity of the quadrupole, this calibration was useful only for 
the quadrupole internally. Since the omegatron calibration is believed to be 
valid, all quadrupole measurements were corrected by a constant factor so that 
the molecular nitrogen density values for the omegatron and the quadrupole, were 
identical in magnitude at 239 km. The internal consistencies of the nitrogen and 
the molecular oxygen data and extensive checking of an identical system in the 
laboratory led to the conclusion that a shift in the sensitivity of the instrument 
occurred that could be represented by a constant. The sensitivities for atomic 
and molecular oxygen were computed with the aid of published ionization cross 
section ratios (Fite and Brachnaarm, 1959) (Kieffer and Bunn, 1966). The accu- 
racy of the cross section ratios is not well known but the errors are probably 
not larger than 1.20 percent. 

No definite explanation could be found ior the apparent sensitivity change 
which seemed to affect the data near minimum angle of attack. The most prob- 
able causes for this orientation-dependent change would be an interference be- 
tween the sun and the ion source, as frequently encountered in TP experiments, 
or possibly an effect of the earth's magnetic field on the gain of the secondary 
electron multiplier. A similar but smaller gain change was observed in a labo- 
ratory test on an identical secondary electron multiplier. However, the data 
obtained at 90 deg angle of attack were in good agreement, both upleg and down- 
leg, and plots of the instrument orientation with respect to the earth’s magnetic 
field vector versus time support this argument. Near the peak, the downleg data 
includes a slight density increase compared with the upleg data. This deviation 
is clearly due to noise and telemetry interference. When the omegatron was 



pointing so as to be illuminated by the sun near apogee the electrometer ampli- 
fier of the omegatron system saturated, driving the telemetry subcarrier oscil- 
lator out of band which generated interference in the adjacent quadrupole channel. 

It was pointed out earlier that little data are available on gas -surface inter- 
actions which could be applied to the analysis of the quadrupole data. Thus , 
errors in the absolute ambient particle densities might occur because of the un- 
certainties in the accommodation coefficients and surface adsoi’ption. Some of 
the uncertainties were eliminated by using data taken at different angles of at- 
tack (minimum and 90 deg), an approach adequate for molecular nitrogen and 
molecular oxygen. The computation using either the constant coefficients a 2 and 
a 3 or the accommodation efficiency e did not yield significantly different results 
because these approaches were alternate ways of accounting for thermal accom- 
modation. In the first attempts to analyze these data specularly reflected par- 
ticles we re included, and the results obtained were only slightly different from 
results obtained without specular reflections. 

The atomic oxygen data required additional critical consideration, for the 
computations wei-e complicated by surface adsorption. It was clear that the as- 
sumed adsorption model did not adequately represent the detailed physical 
processes, e.g., steady state assumption could not fully apply. The final results 
would be similar, however, since the final profile must fall between the extremes 
shown in Figures 13 and 14. Although the data obtained at the times of minimum 
angle of attack allowed the ambient profile to fall into a level which could be un- 
certain by a factor of four, the data obtained at 90 deg angle of attack permitted 
an uncertainty of only approximately 16 percent. This uncertainty did not depend 
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on the surface adsorption model which in our computations narrows the final 
choice only in a somewhat arbitrary way. 

Finally, one might suspect that the change in slope of the ambient density 
profile based on 90 deg angle of attack data might point to a failure of the ion 
source to focus energetic particles. As expected, it is more difficult to focus 
particles into the analyzer section at the 90 deg angle of attack than at small 
angles of attack. However, a careful analysis showed that the focusing correc- 
tion required to bring the ambient density profile, based on 90 deg angle of at- 
tack data, into agreement with the profile based on the minimum angle of attack 
data led to unrealistic assumptions about the electro-dynamic properties of the 
ion source, thus reinforcing our confidence in the focussing properties of the 
ion source. It would also have contradicted the results obtained for molecular 
nitrogen and molecular oxygen where no focusing correction was required but 
where higher particle kinetic energies were encountered. A possible change in 
the focusing of the ion source was therefore discarded. 

The ratios of atomic oxygen density to molecular nitrogen and molecular 
oxygen density are higher than the values measured by some experimenters 
(Hedin, et al., 1964; Hedin and Nier, 1966), but lower than the values measured 
by Mauersberger, et al., 1968, von Zahn, 1970. The preliminary results of the 
recent measurement made with a liquid helium cooled ion source of Offermann 
and von Zahn, 1971 indicate an only slightly higher 0/0 2 ratio when compared 
with extrapolated values of our data. There is, however, relatively good agree- 
ment in the total mass density versus altitude with the corresponding theoretical 
model of the CIRA 1965, COSPAR International Refere nce Atmosphere . For 
comparison the total mass density versus altitude was plotted in Figure 19, 
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together with the value of the CIRA 1965 models three and four. These results 
correspond to a 10.7 cm solar flux of 100 x Iff 22 W/m 2 Hz and 125 x ICC 22 
W/m 2 Hz respectively. The three month average value preceding the launch 
data, August 26, 1966, was 103 x 10 -22 w/m 2 Hz. The geomagnetic index a p was 
6. The model values of the atomic oxygen, molecular oxygen, and molecular 
nitrogen densities were compared in Figure 16 with the measured values and 
good agreement was found between the atomic oxygen and the molecular nitrogen 
densities. The molecular oxygen densities were lower than predicted. Since the 
CIRA model was derived from the total mass density determined from the satel- 
lite drag analysis the close agreement between the present results and the model 
density of molecular nitrogen was striking and it represents an exception from 
results which were obtained for other molecular nitrogen measurements. These 
were consistently lower than the predicted value of the theoretical model. It is 
also interesting to note that the measured ambient temperatures were not in 
agreement with but lower than the model values, a fact which is consistent with 
earlier findings of the omegatron experiment. 
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5. SUMMARY AND CONCLUSIONS 


We report in this paper the results of a different approach to the measure- 
ment of atomic oxygen in the lower thermosphere using a mass spectrometer. 
The primary advantage of the method is the capability of the mass spectrometer 
employed to ionize, focus and analyze a significant number of atmospheric neu- 
tral particles that have not experienced surface collisions and thus are truly 
representative of the ambient conditions. Because of the presence of surfaces, 
the relative velocity between the instrument and the atmosphere and a varying 
angle of attack the measurement will always involve some particles which have 
experienced surface collisions, the number varying strongly with angle of attack 
In this particular application, the ratio of these two ’'classes" of particles varies 
substantially, thus permitting evaluation of the effects of surface collisions. A 
comparison is made between data taken when particles which have experienced 
surface collisions predominate over those which have not, and in circumstances 
where the opposite holds true. As a result it is possible to evaluate the degree 
of thermalization and the degree of absorption of the particles, two properties 
which have strong bearing on the interpretation of data taken during a rocket 
flight. These observations are consistent with the hypothesis that atomic oxy- 
gen is either completely absorbed or, at best, only partially thermalized upon 
impact with instrument surfaces , thus leading to higher concentrations of 0 than 
would be the case if one assumed complete thermalization. 

The results presented here, when compared with data obtained by other 
workers and the assumptions they were obliged to make, seem to verify the 
above assertions. The measured atomic oxygen concentrations are greater, in 
general, than those obtained by the University of Minnesota Group of Prof. A. O. 
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Nier (Hedin et al., 1964, Hedin et aL, 1966, Kasprzak et al. t 1968, and Krankowsky 
et al., 1968) but are less than those measured by the University of Bonn Group of 
Prof. U. von Zahn (Manersbeze et al., 1968, Bitterbeg et al., 1970). At the same 
time, the total density compares favorably with the model atmosphere. It must 
be recognized that all rocket borne atomic and molecular oxygen measurements 
made so far provide only a random sample of a region of the atmosphere in 
which the composition is expected to be considerably more variable than pre- 
dicted by our present atmospheric models. An accurate comparison of the 
various instruments to measure atomic oxygen can therefore only be obtained 
by either a coordinated measurement approach at the lower altitude, for example 
near 120 km, or by designing the instruments so that they reach an altitude where 
atomic oxygen is the major constituent and where the density data can be com- 
pared to satellite drag derived density values. 
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Final Results of the Quadrupole Experiment 
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Figure 3„ Sample of the Telemetry Record Obtained from the Quadrupol 
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Figure 5. Ion source nitrogen density versus flight time for minimum angle of attack, 
90 deg angle of attack, and maximum angle of attack. 
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Figure 6. Ion source molecular oxygen density versus flight time for minimum angle of attack, 
90 deg angle of attack, and maximum angle of attack. 
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Figure 7. Ion source atomic oxygen density versus flight time for minimum angl 
90 deg angle of attack, and maximum angle of attack. 
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Figure 9. Computed molecular nitrogen density versus altitude for minimum angle of attack and 
90 deg angle of attack assuming no thermal accommodation (curves A and B) ( and 
complete thermal accommodation (curves C and D). 


42 




AMBIENT NUMBER DENSITY (CM' 3 ) 


MOLECULAR NITROGEN DENSITY PROFILE 
NASA 18.05 QUADRUPOLE 
WALLOPS ISLAND, VA. 

AUGUST 26, 1966 I4>3I EDT 


io lc M- 



a 2 =o.56, a 3 =o.44 

• UPLEG, BASED ON OUTPUT AT MINIMUM 
ANGLE OF ATTACK. 

+ UPLEG, BASED ON OUTPUT AT 90° 
ANGLE OF ATTACK. 

0 DOWNLEG, MINIMUM ANGLE OF ATTACK. 
0 DOWNLEG, 90° ANGLE OF ATTACK. 

£ =0.52 

□ UPLEG, MINIMUM ANGLE OF ATTACK. 
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Figure 10. Ambient molecular nitrogen density versus altitude obtained from the curves of 
Figure 9 by varying the coefficients aj and a-j or e. The dashed line shows 
the nitrogen density profile obtained from the omegatron. 
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MOLECULAR OXYGEN DENSITY PROFILES 
NASA 18.05 QUADRUPOLE 
WALLOPS ISLAND, VA. 

AUGUST 26, 1966 14 31 EDT 

A : BASED ON PEAK OUTPUT, NO THERMAL ACCOMMODATION 
B= BASED ON OUTPUT AT 90° ANGLE OF ATTACK, NO 
THERMAL ACCOMMODATION 
CEASED ON OUTPUT AT 90° ANGLE OF ATTACK, 
COMPLETE THERMAL ACCOMMODATION 
D ; BASED ON PEAK OUTPUT, 

COMPLETE THERMAL ACCOMMODATION 


ALTITUDE (KM) 


Figure 11. Computed molecular oxygen density versus altitude for minimum angle of attack and 
90 deg angle of attack assuming no thermal accommodation (curves A and B), and 
complete thermal accommodation (curves C and D). 
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MOLECULAR OXYGEN DENSITY PROFILE 
NASA 18.05 QUADRUPOLE 
WALLOPS ISLAND, VA. 

AUGUST 26, 1966 14-31 EOT 


a 2 =o.2i, a 3 =o.79 

♦ UPLEG, BASED ON OUTPUT AT MINIMUM 
ANGLE OF ATTACK. 

+ UPLEG, BASED ON OUTPUT AT 90“ 
ANGLE OF ATTACK. 

o DOWNLEG, MINIMUM ANGLE OF ATTACK. 
(} DOWNLEG, 90° ANGLE OF ATTACK, 

£ - 0.77 

□ UPLEG, MINIMUM ANGLE OF ATTACK 
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Figure 12. Ambient molecular oxygen density versus altitude, obtained from the curves of 
Figure 11 by varying the coefficients 02 and or 6. 
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re 13. Computed atomic oxygen density versus altitude for minimum angle of attack assumi 
(1) complete adsorption of particles striking the surface (curve A), (2) no adsorption 
and no thermal accommodation of particles striking the top deck (curve B), (3) no 
adsorption and complete thermal accommodation of particles striking the top deck 
(curve C). A theoretical profile computed by assuming diffusive and thermal 
equilibrium and by using the omegatron temperature is shown by the dashed line. 
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Figure 14. Computed atomic oxygen density versus altitude for 90 deg angle of attack assuming 
(1) complete adsorption of particles striking the surface (curve A), (2) no adsorption 
and no thermal accommodation of particles striking the top deck (curve B), (3) no 
adsorption and no thermal accommodation of particles striking the top deck (curve C). 
The two dashed lines each result from assuming diffusive and thermal equilibrium 
and from using the omegatron temperature. The top dashed line assumes complete 
adsorption at 170 km. The bottom dashed line assumes complete accommodation 
and no adsorption at high altitudes. 
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Figure 15. Ambient Atomic Oxygen Density Versus Altitude 
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Figure 16. Ambient particle densities of molecular nitrogen, atomic and molecular 
oxygen versus altitude. 
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Figure 18. Ambient temperature versus altitude. The dashed lines show the value of the CIRA 
1965 models 3 and 4 for 1400 hours. 


51 



MASS DENSITY (KG/M 3 ) 






APPENDIX A 


A SIMPLE MODEL FOR THERMAL ACCOMMODATION OF 
NEUTRAL PARTICLES ON THE SURFACE OF THE ION SOURCE 


The derivation of the thermal accommodation equation is based on the as- 
sumption that all particles which collide with the top surface of the ion source 
lose an amount of energy which is proportional to the initial energy with which 
they strike the surface. 

Defining an accommodation efficiency as 
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: 1 1 2 


C 2 

j 2 3 

C 2 

s 


u-v 


u-v 


(A.l) 


where 

u = probe velocity, 

u-v - initial velocity before the collision, 
v' = final velocity after the collision, 

C 0 = j/2kT fi /in, most probable speed in a gas at ambient temperature, 
C s = y'2kT s /m, most probable speed in a gas at the temperature of the 
surface, 

where 

k = Boltzmann's constant, 1.38 x 10 -23 J/°K, 

T 0 = ambient neutral particle temperature, 

T s = surface temperature, 
m = molecular weight of the specie measured. 


A-l 



From Equation (A.l), we obtain the final velocity ratio 
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(A .2) 


The Maxwell-Boltzmann velocity distribution which describes the ambient gas in 
a moving coordinate system is 


f(u,v) - exp - 

V C 0 3 7T3/2 


U- V 


2 



(A. 3) 


With the assumption made above, the term |u~vj 2 /C 0 2 in the exponent of Equa- 
tion (A. 3) is replaced by Equation (A. 2) and the distribution function of the out- 
going gas can be written 


V’ 
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|u -v | 2 

|v| 2 " 

s 

= K exp - 

(1 - e) c* 

o 

+ £ 

C 2 

s 


f* (u, v) - K exp - 

where K is a constant to be determined from the condition 

jf 1 (u,v) d 3 v = 1. 

(all velocity space) 

This condition yields 
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C 3/2 77-3/2 


3/2 


1 + T " 1 £ 


and, therefore, the velocity distribution of the outflowing gas is 
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Equating r t to r a , we have 


n o C 0 


f(S z ) 


n s C 0 




and, solving for n s , 
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f(S z ) 

n o H(€,S z ) • 


(A.8) 


The density in the ionizing region is now obtained by dividing the outgoing 
differential flux by the particle velocity, by integrating over all possible veloci- 
ties, and then by multiplying by the solid angle subtended by the surfaces. If n d 
is the particle density in the ionization region, 


n, 


fi s 
27 r. 


/ dl 2 ~ao pao ~oo 

= f, (n,v) dv x dv y dv 2 

z Jq J- co J-CO 


(A. 9) 


where Q s is the solid angle subtended by the surface, 

Q s - 27 t(cos 0 2 - cos0 x ) . 


Thus, after integration, 

(cos 6 ^ ~ cos 6 j) 

n d = 2 n s [ 1 + erf ( S z)] ex P ( _S i z ) ( A ' 10 ) 


and by substituting Equations (A. 7) and (A.8) into (A. 10), 
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(A.4) 


Proceeding as in Niemann and Kreick (1966) by computing and equating the in- 
coming and outgoing flux we find the surface particle density. 

The incoming flux is given as 


r, 


n. 


*•00 s%CX) /-• CO 

0 J — co * 2 - 


n o C 0 


V z f(u,v) dv x dv y dv z - f (S z ) , (A. 5) 


where 

n 0 = ambient particle density, 

C 0 - average thermal velocity, C Q = ]/skT 0 firm , 
f (S) = exp (-S 2 ) + W S[ l + erf (S)] , 

S z = u z /C Q , the ratio of the z-component of the vehicle velocity to the 
most probable thermal velocity. 

Similarly, the outgoing flux is 
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J n J -CO J - 
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V z f' (u.v) dv x dv y dv z - -J- H(e,s z ) , (A. 6) 


where n s is the surface gas density and 
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with 
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APPENDIX B 


THE ATOMIC OXYGEN-SURFACE INTERACTION MODEL 

Atomic oxygen was partially adsorbed on the deck of the quadrupole. The 
general features of this adsorption were investigated by using the simple model 
derived below. This model assumes that (1) the time required to achieve steady 
state is short compared to the tumble period, (2) the model parameters are 
functions of the quadrupole angle of attack, that is, functions of the angular dis- 
tribution of the incoming flux, and (3) the second and all higher adsorbed gas 
layers are quite similar to ea.ch other but different from the first layer. 

Let 

b t = evaporation rate from the first layer, 
b = evaporation rate from the second and higher order layers, 
a t - sticking probability on unoccupied sites, 
a = sticking probability on occupied sites, 
y l = recombination coefficient of atomic oxygen for the first layer 
y - recombination coefficient of atomic oxygen for the second and higher 
order layers, 

c j = rate of diffusion of atomic oxygen into the solid, 

Fj = incoming atomic oxygen flux, 

P 2 = outgoing atomic oxygen flux, 

S i = fractional surface coverage of the i th surface containing i layers of 
atomic oxygen, i - 0, 1, . . . 

Assuming steady state surface conditions, we obtain from the conservation of 
fractional coverage for the first layer, 


ZB- 


B-l 



S i b i + r i 7i s i + c i( 1 " S o) ~ r i a i s o- 


(B.l) 


Similarly for the second layer, 

r i a i S o + S 2 b + y s 2 = bj + r J, y 1 s x + Cj(l - s 0 ) + Fj a s t . 

From Equations (B.l) and (B.2), 

S 2 b + r i 7 S 2 - Tj a Sj = 0 . 

Thus, in general we can write for the second and higher order layers 
b + 7 S i - Pj a S i _ 1 - 0, i = 2, 3, . . . . 


Solving Equation (B.l) for S 3 . 
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and solving Equation (B.4) for S i we have 
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where 
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(B. 6 ) 



X = 


7 r 1 + b ’ 

From the definition of S i we have the condition 

f S i = 1, (B.7) 

i = 0 

or with Equations (B.5) and (B.6) 

s 0 + Si + s 2 + . . . 

- S 0 + (YS 0 - Z) + (YS 0 - Z) X + (YS 0 - Z) X 2 + . . . 


- S Q + (YS 0 - Z) £ X* = 1 . 

i =0 


(B.8) 


If X <1, then the series converges and 


1 - X + Z 
S o - 1 - X + Y 


(B.9) 


The above equation is valid if the surface properties of atomic oxygen are 
invariant from layer to layer. Now, 


r 2 - r x - - r l7 s,[ - Cl (i - s 0 ) , 


and if the fractional adsorption is defined 
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r - r 
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= 7 ! Si + 7(1 - S 0 - s x ) + — (1 - S 0 ) , (B.10) 


after rearranging 
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by using Equation (B.5) for S 1 , (B.9) for S 0 , and by rearranging, 
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(B.12) 


Equation (B.12) reduces to 



( Cl \l 
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(y, - y) (1 - X) + (y + r j 


(1 - X + Y) 


(B.13) 


Substituting the values of X, Y, and Z into Equation (B.13) yields the general 
form of the adsorption equation of the assumed model: 
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The coefficients of Equation (B.14) are obtained by fitting the data to a 
theoretical profile so that they are consistent with the assumptions made in the 
adsorption model. The simplest consistent form is obtained by letting c 1 = 
y, = 0. With this simplification, 
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By using the above assumptions, for the minimum angle of attack, 
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and for the 90 deg angle of attack, 
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FIGURE CAPTIONS 


Figure 1. The thermosphere probe in the clam shell nose cone. 

Figure 2. Cross section of the quadrupole ion source and lens chamber with 
the electron beam trajectory. 

Figure 3. Sample of the telemetry record obtained from the quadrupole. 

Figure 4. (a) Maximum and minimum angle of attack between the thermosphere 
probe axis and the velocity vector versus time from launch. 

(b) Expanded view of a section of Figure 4A. 

Figure 5. Ion source nitrogen density versus flight time for minimum angle of 
attack, 90 deg angle of attack, and maximum angle of attack. 

Figure 6. Ion source molecular oxygen density versus flight time for minimum 
angle of attack, 90 deg angle of attack, and maximum angle of attack. 

Figure 7. Ion source atomic oxygen density versus flight time for minimum 
angle of attack, 90 deg angle of attack, and maximum angle of attack. 

Figure 8. Simplified model of the geometry of the quadrupole ion source. 

Figure 9. Computed molecular nitrogen density versus altitude for minimum 
angle of attack and 90 deg angle of attack assuming no thermal accommo- 
dation (curves A and B), and complete thermal accommodation (curves C 
and D). 



Figure 10. Ambient molecular nitrogen density versus altitude obtained from 
the curves of Figure 9 by varying the coefficients a 2 and a 3 or e. The 
dashed line shows the nitrogen density profile obtained from the omegatron. 

Figure 11. Computed molecular oxygen density versus altitude for minimum 
angle of attack and 90 deg angle of attack assuming no thermal accommoda- 
tion (curves A and B), and complete thermal accommodation (curves C 
and D). 

Figure 12. Ambient molecular oxygen density versus altitude, obtained from the 
curves of Figure 11 by varying the coefficients a 2 and a 3 or e. 

Figure 13. Computed atomic oxygen density versus altitude for minimum angle 
of attack assuming (1) complete adsorption of particles striking the surface 
(curve A), (2) no adsorption and no thermal accommodation of particles 
striking the top deck (curve B), (3) no adsorption and complete thermal ac- 
commodation of particles striking the top deck (curve C). A theoretical 
profile computed by assuming diffusive and thermal equilibrium and by 
using the omegatron temperature is shown by the dashed line. 

Figure 14. Computed atomic oxygen density versus altitude for 90 deg angle of 
attack assuming (1) complete adsorption of particles striking the surface 
(curve A), (2) no adsorption and no thermal accommodation of particles 
striking the top deck (curve B), (3) no adsorption and no thermal accommo- 
dation of particles striking the top deck (curve C). The two dashed lines 
each result from assuming diffusive and thermal equilibrium and from using 
the omegatron temperature. The top dashed line assumes complete 



adsorption at 170 km. The bottom dashed line assumes complete accommo- 
dation and no adsorption at high altitudes. 

Figure 15. Ambient atomic oxygen density versus altitude. 

Figure 16. Ambient particle densities of molecular nitrogen, atomic and 
molecular oxygen versus altitude. 

Figure 17. Ambient particle density ratios of atomic to molecular oxygen, 

atomic oxygen to molecular nitrogen, and of molecular oxygen to molecular 
nitrogen versus altitude. 

Figure 18. Ambient temperature versus altitude. The dashed lines show the 
value of the CIRA 1965 models 3 and 4 for 1400 hours. 

Figure 19. Total mass density of major constituents versus altitude. The 

dashed lines show the value of the CIRA 1965 models 3 and 4 for 1400 hours. 

Table I. Final results of the quadrupole experiment showing number densities 
(n a (N 2 ), n a (0 2 ), n a (0)), number density ratios (0/0 2 , 0/N 2 , 0 2 /N 2 ), mass 
density of major constituents (p), and mean molecular weight of major 
constituents (M). 




